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Abstract. Exact models of uniformly rotating strange 
stars, built of self bound quark matter, are calculated 
within the framework of general relativity. This is made 
possible thanks to a new numerical technique capable to 
handle the strong density discontinuity at the surface of 
these stars. Numerical calculations are done for a simple 
MIT bag model equation of state of strange quark mat- 
ter. Evolutionary sequences of models of rotating strange 
stars at constant baryon mass are calculated. Maximally 
rotating configurations of strange stars are determined, as- 
suming that the rotation frequency is limited by the mass 
shedding and the secular instability with respect to ax- 
isymmetric perturbarions. Exact formulae which give the 
dependence of the maximum rotation frequency, and of the 
maximum mass and corresponding radius of rotating con- 
figurations, on the value of the bag constant, are obtained. 
The values of T/W for rapidly rotating massive strange 
stars are significantly higher than those for ordinary neu- 
tron stars. This might indicate particular susceptibility of 
rapidly rotating strange stars to triaxial instabilities. 
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1. Introduction 

A deconfined, beta stable quark matter is composed of the 
u, d, and s quarks, and in contrast to "ordinary" baryon 
matter has the strangeness per unit baryon number ~ — 1 
(therefore the name 'strange quark matter'). An intriguing 
possibility that such strange matter could be the absolute 
ground state of matter at zero pressure and temparature 
was pointed out in an influential paper of Witten (1984) 
[such a possibility has been contemplated before by Bod- 
mer (1971)]. This possibility is not excluded by what we 
know from laboratory nuclear physics. If the hypothesis of 
strange matter is true, then some of neutron stars could 
be actually strange stars, built entirely of strange matter 
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(Haensel et al. 1986, Alcock et al. 1986; arguments against 
the existence of strange stars have been presented by Al- 
par (1987), and Caldwell & Friedman (1991) ). 

Detailed models of nonrotating strange stars were con- 
structed in (Haensel et al. 1986, Alcock et al. 1986). Cal- 
culations concerning physics and astrophysics of strange 
stars, published up to 1991, were reviewed in (Madsen & 
Haensel 1991). References to more recent work on strange 
stars can be found in review articles by Weber et al. (1995) 
and Madsen (1999), and in the monograph of Glcnden- 
ning (1997). Most recently, strange stars have been in- 
voked in models of X-ray bursters (Bombaci 1997, Cheng 
et al. 1998) and of gamma-ray bursters (Dai & Lu 1998). 

Rapid rotation of strange stars has become a topic 
of interest in 1989, after a sensational (but subsequently 
withdrawn) claim of detection of a 0.5 ms pulsar in 
SN 1987A. In several papers, strange star has been ad- 
vanced as an appropriate model for a half-millisecond 
pulsar (Glendenning 1989a, b, Frieman & Olinto 1989). 
However, several authors pointed out that the experi- 
mental constraints on the strange matter model imply 
that strange stars could not rotate so fast (Haensel and 
Zdunik 1989, Lattimer et al. 1990, Zdunik & Haensel 1990, 
Prakash et al. 1990). Subsequent studies of rapid rota- 
tion of strange stars focused on the effect of rotation on 
the structure of the normal crust (Glendenning & Weber 
1992). 

Recently, it has been shown that fast rotation of hot, 
young neutron stars is severely limited by the emission of 
gravitational radiation, due to r-mode instability [for a re- 
view, see e.g. Stergioulas (1998)]. However, in contrast to 
newly born neutron stars, hot strange stars are not sub- 
ject to the r-mode instability (Madsen 1998). Therefore, 
strange stars formed in collapse of rotating stellar cores 
can rotate very fast. This gives additional motivation for 
studying rapid rotation of strange stars. 

It should be stressed, that nearly all calculations of uni- 
formly rotating models of strange stars were done within 
the slow rotation approximation (see, e.g., Colpi & Miller 
1992, Glendenning & Weber 1992). Actually, in the calcu- 
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lations of Glendenning & Weber (1992) the slow rotation 
scheme, pioneered by Hartle (1967, 1973) and Hartle & 
Thorne (1968), was supplemented by a "self-consistency 
condition" at the Keplerian frequency (Weber & Glen- 
denning 1991, 1992). However, the slow rotation approx- 
imation is not valid near the shedding (Keplerian) limit. 
As discussed by Salgado et al. (1994a, b), the improved, 
"self-consistent" version of this approximation, overesti- 
mates by more than 10% the maximum rotation frequency 
of neutron star models. 

So far, the only exactf] calculations of the rapidly rotat- 
ing models of strange stars were done by Friedman (quoted 
in Glendenning 1989a, b) and Lattimer et al. (1990). How- 
ever, the precision of their numerical results, calculated 
for one specific equation of state of strange matter, was 
not checked using an independent exact calculation. Also, 
the fact that their 2-D code was based on the Butterworth 
& Ipser (1976) scheme, rises some suspicion concerning its 
precision when applied to the equation of state of strange 
matter, characterized by huge density discontinuity at the 
stellar surface. Indeed Butterworth & Ipser (1976) con- 
sidered the case of a density discontinuity at the stellar 
surface - they computed configurations of rotat ing homo- 
geneous bodies, but as we discuss below (Sec. 4.3), they 
treated properly the discontinuity only in the radial di- 
rection, which is not sufficient for objects which deviate 
substantially from the spherical symmetry. 

In the present paper we perform exact calculations 
of rapidly rotating strange stars, composed entirely of 
strange quark matter. Our calculations are based on 
a multi-domain spectral method recently developed by 
Bonazzola et al. (1998b). Such a multi-domain technique 
enables us to treat exactly the density discontinuity at the 
surface of strange stars. The calculations are performed for 
a family of equations of state of strange quark matter. On 
one hand, precision of our calculation is checked using in- 
ternal error estimators. On the other hand, the validity of 
scaling relations for the rotating configurations, displayed 
by our numerical results, yields an additional test of pre- 
cision of our calculations. Finally, we address the problem 
of triaxial instabilities of rapidly rotating strange stars, 
and point out differences with respect to normal neutron 
stars. 

In Sect. H we discuss the equation of state of strange 
matter, based on the MIT bag model of quark matter. The 
set of equations to be solved and the numerical procedure 
is presented in Sect. ^. Stationary, uniformly rotating con- 
figuration of strange stars are studied in Sect. |J, where we 
also derive exact scaling relations for the parameters of 
the rotating strange star models. Parameters of the max- 
imally rotating configurations, stable with respect to the 
axially-symmetric perturbations, are calculated in Sect. B, 



and an exact formula for the maximum rotation frequency 
of rotating strange stars is derived. We derive also exact 
formulae for mass and radius of the maximum mass config- 
uration of rotating strange stars. In Sect. |[ we study the 
problem of triaxial instabilities of rapidly rotating strange 
stars, and show, that rapidly rotating strange stars might 
be more susceptible to these instabilities than ordinary 
neutron stars. Finally, Sect, ^contains a discussion of our 
results and the conclusion. 



2. Equation of state and static models of bare 
strange stars 

Equation of state of strange matter will be based on the 
MIT bag model. Baryon number density of strange matter 
is n — -|(n u + + n s ), where n u is the number density of 
the u-quarks etc. In what follows, we will use the simplest 
model of self-bound strange quark matter, ignoring the 
strange quark mass and neglecting the quark interactions 
except for the confinement effects described by the bag 
constant. We consider thus massless, noninteracting u, d, 
s quarks, confined to the bag volume. In the case of a 
bare strange star, the boundary of the bag coincides with 
stellar surface. 

Using the model described above, one can express the 
energy density, e, and the pressure, P, of strange quark 
matter, as functions of the baryon number density, n, in 
the following form 



e = an 3 + B , 

1 i 

P = -an 3 — B 
3 



(1) 



where B is the MIT bag constant, and the parameter a is 
given by 



a = -nihe = 952.371 MeV fm 
4 



(2) 



The bag constant B describes the difference in the energy 
density of the true (real) vacuum and that of the QCD vac- 
uum. This parameter plays a crucial role in the MIT bag 
model, being actually responsible for the quark confine- 
ment. Let us mention, that in the case of a bare strange 
star of ~ M Q we are dealing with a bag of a radius of 
~ 10 km, containing ~ 10 57 quarks. 

The constant B represents a natural unit for both the 
energy density and the pressure. Defining dimcnsionlcss 
energy density and pressure of strange matter, e = e/B, 
P = P/B, we can write down a dimensionless form of the 
EOS of strange matter, 



P=^-4). 



(3) 



1 In this article, the term "exact" is relative to the treatment 
of rotation and is used to distinguish from the slow rotation 
approximation 



The above dimensionless form of the EOS of strange mat- 
ter constitutes a basis for the derivation of the scaling laws 
relating the families of static models of strange stars, cal- 
culated for different values of B (Witten 1984, Haensel et 
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Fig. 1. Adiabatic index of strange quark matter, 7, versus 
baryon density, n. Two vertical dashed lines correspond 
to strange star surface (left line), and maximum density 
within static strange star models (right line). Horizontal 
dashed line corresponds to 7 = 4/3, characteristic of free 
ultrarelativistic Fermi gas. 

al. 1986). For the (gravitational) mass, baryon mass, and 
radius of static configuration with maximum allowable 
mass we get (assuming G = 6.6726 10~ n m 3 kg _1 s~ 2 , 



M Q = 1.989 1 30 kg and c : 
1.9638 



2.99792458 10 8 ms" 1 ): 



71 stat 
iVi B.max 
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2.6252 
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! 



km 



(4) 



3 ). The baryon mass is de- 
Nb, where Ab is the total 



where B 60 = 5/(60 MeV fm 
fined as M B = 1.66 10~ 27 kg 
baryon number of stellar configuration. 

Consider the case of strange matter at zero pressure 
(such a situation corresponds to the surface of a "bare 
strange star"). For massless, noninteracting quarks, we 
obtain the following set of parameters, characterizing the 
properties of strange matter at zero pressure, 



n = 0.28665 (S 60 ) 3/4 fm~ 3 , 

po = e /c 2 = 4.2785 x 10 17 B m kg ixT 3 



(5) 



Strange matter would be the real ground state of matter 
at zero pressure, only if the energy per unit baryon num- 
ber, at zero pressure, Eq = e /n = 837.26(£?6o) 1 ^ 4 MeV, 
was below the energy per baryon in the maximally bound 
terrestrial nucleus, 56 Fe, equal 930.4 MeV. This implies 



£?6o < 1-525. On the other hand, a consistent model 
should not lead to a spontaneous fusion of neutrons into 
droplets of u, d quarks (which would eventually trans- 
form into droplets of strange matter - "strangelets"). In 
view of the fact, that a uc i/a = (1 + 2 4 / 3 )/3, we get 
E ( ud] = 1.127^ uds) . Therefore, the condition of the sta- 
bility of neutrons with respect to a spontaneous fusion 
into strangelets E ( Q ud) = 988.5( J B 60 ) 1/4 MeV > 939.6 MeV, 
which implies Bqq > 0.9821. Summarizing, within the sim- 
plest model of strange matter, the bag constant is con- 
strained by 0.9821 < B m < 1.525. 

An important quantity, relevant for the pulsations and 
stability of strange stars is the adiabatic index of strange 
matter, defined as 7 = (n/P)dP/dn. Dependence of 7 on 
baryon density of strange matter is shown in Fig. |l|; it is 
qualitatively different from 7(71) for ordinary neutron star 
matter. The values of 7 in the outer layers of strange stars 
are very large. Even at the highest densities, allowed for 
strange stars, the value of 7 is significantly higher, than 
the ultrarelativistic Fermi gas value 4/3, predicted by the 
asymptotic freedom of QCD at n =>• 00. 

3. Formulation of the problem and numerical 
code 

3.1. Equations of stationary motion 

We refer to Bonazzola et al. (1993) for a complete descrip- 
tion of the set of equations to be solved to get general rela- 
tivistic models of stationary rotating bodies. Let us simply 
recall here that, under the hypothesis of stationarity, ax- 
isymmetry and purely azimuthal motion (no convection), 
a coordinate system (t, r, 8, (p) can be chosen so that the 
spacetime metric takes the form 

ds 2 = -A 2 dt 2 + B 2 r 2 sin 2 6{dip ~ N v dt) 2 



+A 2 (dr 2 +r 2 dd 2 ) , 



(6) 



where A, N v , A and B are four functions of (r,8). The 
Einstein equations result in a set of four elliptic equations 
for these metric potentials^]: 

4:nA 2 (E + 3P+(E + P)U 2 ) 



+ B 2 r 2 sm 2 d {dNV)2 _ ^ + ^ 



2A 2 
A 3 (A^rsinfl) 



NA 2 . 

-167T (E- 

B v 



P)U 



~r sin 6 dN v 9(3/3 — v) 
A 2 [(NB - 1) r sin 6] = 16nNA 2 BPr sin 9 
A 2 (v + a) = 8ttA 2 [P+ [E + P)U 2 ] 
3B 2 r 2 sin 2 9 



AN 2 



■(dN*) 2 - {dvf 



(7) 

(8) 
(9) 

(10) 



where the following abreviations have been introduced: 
v:=\nN , a:=\nA, (3:=\nB, (11) 



in this section G = c 
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(12) 
(13) 
(14) 
(15) 



In Eqs. (Q)-(JT^), U and E are respectively the fluid 3- 
velocity and energy density, both as measured by the lo- 
cally non-rotating observer: U — Br sin 9(Cl — N ip )/N, 
E = T 2 (e + P)-P,T=(l-U 2 )~ 1 / 2 . 

These equations are supplemented by the first integral 
of motion 



const. 



(16) 



where H is pseudo-enthalpy (or log-enthalpy) defined by 
fP N 



H = \n 



nE Q 



(17) 



The dimensionless EOS of strange matter can be 
parametrized in terms of the H variable as 



e = 3e 4H + 1 
P = c AH - 1 . 



(18) 



Consequently, the baryon number density can expressed 
in terms of H by 



n e 



3H 



(19) 



where for our model the baryon number density of strange 
matter at zero pressure no = (3£?/et) 3 / 4 . 

The equations of stationary motion involve the quan- 
tities with dimensions of energy density and powers of 
length, and universal constants G and c. These equations 
can be rewritten in a dimensionless form, if we introduce 
the dimensionless quantities e, P, r, f2, via following rela- 
tions, 



e 

P 



= Be, 
= BP 



n = 




(20) 



Notice, that the solutions of dimensionless equations of 
stationary motion will not depend explicitly on B. In order 
to recover the solution in conventional units, for a specific 
value of B, one has just to use relations (pfj|). 

3. 2. Numerical procedure 

The non- linear elliptic equations (|7j)-(|io|) are solved it- 
eratively by means of the multi-domain spectral method 



developed recently by Bonazzola et al. (1998b). In this 
method, the whole space is divided into three domains: 
Dl: the interior of the star, D2: an intermediate domain 
whose inner boundary is the surface of the star and outer 
boundary a sphere located at r := R23 ~ 2r oq (where r eq is 
the equatorial coordinate-radius of the star), and D3: the 
external domain whose inner boundary is the outer bound- 
ary of D2 and which extends up to infinity, thanks to the 
compactification u = l/r. A mapping (£, 6') 1— > (r,0) sat- 
isfying 9 = 9' is introduced in each domain so that the 
domain boundaries lie at a constant value of the coordi- 
nate £ (typically, £ = —1, or 1). Explicitly, this mapping 
reads 

inDl: r = ai [£ + (3£ 4 - 2£ 6 ) Fx(6")] , £ e [0, 1] , (21) 



in D2: 



in D3: 



a 2 

1 

r 



£ 3 -3£ 



-F 2 (9') + R23 , 
£e[-l,l], (22) 



2R 



2:S 



-(1-0. £e[-i.i] 



(23) 



where the functions F\{9') and ^2(6*') are related to the 
equation r — S(9) of the stellar surface by S(9) =_ai[l + 
Fi(9)} = a 2 [-2 



F 2 {9)} + i? 23 - The mapping @- 
is a specialization to the axisymmctric case of the 3-D 
mapping introduced in Bonazzola et al. (1998b), to which 
the reader is referred for more details. 

A spectral expansion of each relevant field is then per- 
formed with respect to the coordinates (£, 6'). We use 
Chebyshev polynomials in £ and trigonometric polynomi- 
als or Legendre polynomials in 9' (see Bonazzola et al. 
1999 for the use of spectral methods in relativistic astro- 
physics). 

Since the discontinuities in the physical fields (strong 
discontinuity in the density, cf. Fig. || below, discontinu- 
ity in the second derivative for the metric potentials) are 
located at the boundary between two domains (namely 
Dl and D2), the applied spectral method is free from any 
Gibbs phenomenon and leads to a very high precision. 
This is illustrated in Fig. 5 of Bonazzola et al. (1998b) 
which concerns the case of a rotating constant-density 
body (hence with a strong discontinuity of the density 
field at the surface) in the Newtonian theory. An analyti- 
cal solution is available in this case (MacLaurin ellipsoid) 
and can be used to evaluate the accuracy of the code. The 
code gives a rapidly rotating MacLaurin ellipsoid with a 
relative error of the order 10~ 12 with 49 (resp. 25) de- 
grees of freedom in r (resp. 9). Figure 5 of Bonazzola et 
al. (1998b) also demonstrates that the error is evanescent, 
i.e. that it decreases exponentially with the number of de- 
grees of freedom (i.e. grid points). 

This represents a major improvement with respect to 
the spectral method developed previously by Bonazzola 
et al. (1993) to compute stationary configurations of ro- 
tating bodies in general relativity. This method has been 
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Fig. 2. Meridional plane cross section of a rapidly rotat- 
ing strange star, for B = 60 MeV fm -3 . The star's baryon 
mass is Mb = 2.00 Mq, gravitational mass M — 1.60 M©, 
and rotation period P — 0.87 ms. The coordinates x and z 
are defined from the coordinates (r, 9, (p) introduced in the 
line element ^ by x — r sin 9 cos ip and z — rcos9. The 
rotation axis is the z-axis. The various lines are isocon- 
tours of the log-enthalpy H . Solid lines indicate a positive 
value of H and dashed lines a negative value. In this latter 
case (vacuum), H is defined from Eq. (|l6|). The thick solid 
line denotes the stellar surface. 



used to compute rapidly rotating neutron star models by 
Salgado et al. (1994a, b), Haensel et al. (1995), Goussard 
et al. (1997, 1998) and Nozawa et al. (1998) but it could 
not have been employed as such to compute models of ro- 
tating strange stars: the Gibbs phenomenon induced by 
the density discontinuity at the stellar surface would have 
been too large. 

3.3. Description of a typical model 

A rapidly rotating strange star model obtained from 
the code is shown in Figs. || and [| Assuming B = 
60 MeVfm~ 3 , this model has a baryon mass Mb = 
2.00 M Q , a gravitational mass M = 1.60 Mq, and a ro- 
tation period P = 0.87 ms. The central baryon density is 
n c = 0.42 fm -3 , the central values of the metric potentials 
N, A, B and [cf. Eq. (§)] are respectively 0.65, 1.47, 
1.47 and 0.46 fl. The boundary between the computational 
domains Dl and D2 introduced in Sect. 3.2 coincides with 
the stellar surface (thick solid line in Fig. |2|). The bound- 
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Fig. 3. Profiles of the proper energy density e in the equa- 
torial plane (9 — tt/2) and along the rotation axis (9 = 0) 
for the rotating strange star model shown in Fig. 0. 



ary between the computational domains D2 and D3 is out 
of the scope of the Figure. The energy density profile corre- 
sponding to this model is shown in Fig. ||[ Note the strong 
discontinuity at the stellar surface: the jump in density is 
~ 2/3 of the central density. 

3.4- Treatment of configurations at the Keplerian limit 

The maximum rotation rate that a star can sustain is 
reached when the velocity at the equator equals the ve- 
locity of an orbiting particle (a higher velocity would re- 
sult in a centrifugal break-up of the star). This limit is 
called the Keplerian limit. It is reached when the deriva- 
tive dH/dr vanishes at the stellar equator. The surface 
of the star is then no longer smooth but exhibits a cusp 
along the equator (see Fig. ||). Such a non-diffcrcntiable 
surface cannot be decribed by the mapping (|2^)-(|2^) be- 
cause the functions F\ (9) and Fi (9) are assumed to be ex- 
pandable in cos(k9) series (cf. Bonazzola et al. (1998b)), 
which implies that they are smooth functions of 6. The 
solution to this problem consists in freezing the adapta- 
tion of the mapping to the stellar surface when the ratio 
(dH/dr) cq J (dH / dr) po \ c passes below a certain threshold 
during the iteration process. For instance, this threshold 
was chosen to be 0.28 in the computation shown in Fig. [ij. 
Consequently, in the final result, the density discontinuity 
at the stellar surface no longer coincides with the bound- 
ary between the computational domains Dl and D2, ex- 
cept at the equator. In this case, a Gibbs phenomenon is 
present. The accuracy of the calculation is then lower than 
when the mapping is adapted to the surface of the star. 
However, since the stellar interior covers most of domain 
Dl (cf. Fig. |[), the Gibbs phenomenon is rather limited. 
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Fig. 4. Same as Fig. || but for a strange star rotat- 
ing at the Keplerian limit. The star's baryon mass is 
Mb = 2.34 M Q , gravitational mass M = 1.89 M Q , and 
rotation period P = 0.82 ms. Thin dashed and solid lines 
are isocontours of the log-enthalpy H . Thin dashed lines 
correspond to negative values, and thin solid lines to pos- 
itive values of H, respectively. The thick dot-dashed line 
denotes the boundary between the computational domains 
Dl and D2. 



It is notably less severe than if Dl would have been a 
sphere of radius r oq , as in our previous numerical method 
(Bonazzola et al. 1993). 

3. 5. Tests of the numerical code 

Various tests have been passed by the code. First of all, 
we recover previous results, e.g. those presented in Nozawa 
et al. (1998), when using an EOS for neutron star mat- 
ter instead of the strange quark matter EOS presented in 
Sect. H 

Regarding the treatment of the density discontinuity at 
the stellar surface, the code has been tested on the com- 
putation of MacLaurin ellipsoids (homogeneous rotating 
bodies in the Newtonian regime), and gives excellent re- 
sults, as recalled in Sect. 3.2: the relative accuracy reaches 
1(T 12 ! 



Another type of test is the evaluation of the virial iden- 
tities GRV2 (Bonazzola 1973, Bonazzola & Gourgoulhon 
1994) and GRV3 (Gourgoulhon & Bonazzola 1994), this 
latter being a relativistic generalization of the classical 
virial theorem. GRV2 and GRV3 are integral identities 



which must be satisfied by any solution of the Einstein 
equations (|7|)-(fL0|) and which are not imposed during the 
numerical procedure (cf. Nozawa et al. 1998 for details on 
the computation of GRV2 and GRV3). When presenting 
numerical results, we will systematically give the accuracy 
by which the numerical solution satisfies these virial iden- 
tities. 



As discussed in Sect. 3.4, the computation of configu- 
rations rotating at the Keplerian limit is not free of Gibbs 
phenomenon. In order to gauge the resulting numerical er- 
ror, we performed various computations of the same Kep- 
lerian configuration, by varying the numbers of coefficients 
in the spectral method and by varying the threshold on 
(dH/dr) cq I (dH/dr) po \ c for freezing the mapping. We 
systematically used three sets of numbers of coefficients (= 
number of collocation points): (N r , Ng) = (33, 17), (37, 19) 
and (41,21) in each of the three domains. The threshold 
on (dH/dr) cq / {dH/dr) po \ c was varied from 0.28 down 
to 0.20. The variation of all these parameters lead to a 
relative change of the numerical solution of the order of 
or below 1%. This gives an estimation of the error of our 
method in computing strange stars at the Keplerian limit. 
Note the GRV2 and GRV3 errors for Keplerian configu- 
rations revealed to be better than this. For instance, for 
the configuration shown in Fig. ||, the GRV2 (resp. GRV3) 
error is 3 10~ 5 (resp. 1.3 10~ 3 ). 



4. Numerical results 

4-.1. General properties 

For rigid rotation, the equilibrium configurations span a 
two-dimensional domain. This domain is shown in Fig. 
in the central density - gravitational mass plane. Selected 
configurations are listed in Table [j] The symbols have the 
following meaning: M is the gravitational mass; Mb is the 
baryon mass; Q is the angular velocity, P the correspond- 
ing rotation period; n c is the central baryon density; p c 
is the central proper energy density divided by c 2 ; H c is 
the central log-enthalpy [Eq. (|l7|)] ; i? eq is the circumferen- 
tial radius, i.e. the length of the equator [as given by the 
metric (^|)] divided by 27r; r cq is the equator r coordinate; 
rp ie/Veq is the coordinate oblateness of the star; / is the 
moment of inertia (defined as J/fl); J is the star angular 
momentum; Q is the quadrupole moment defined as in Sal- 
gado et al. (1994a) or Laarakkers & Poisson (1999); T/W 
is the "kinetic to gravitational energy" ratio (see Sect. ||); 
U eq is the rotation velocity at the equator as measured 
by a locally non-rotating observer (those observer whose 
4- velocity is the normal to the t = const hypersurfaces); 
z eq is the redshift for an emission at the equator and in 
the direction of rotation, z^ q the redshift for an emission 
at the equator and in the direction opposite to rotation, 
Zp O i the redshift at the stellar pole; N c , N^, A c and B c 
are the metric potentials [Eq. (0)] at the stellar center. Fi- 
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Table 1. Selected strange star models (B = 60 MeVfm 3 ); symbols are denned in Sect. 



4.1. 



Model 


M = 1.4 Mq 

D. = 


MSS 


M = 1.4 Mq 

P = 1.56 ms 


M rot 

■"'max 


^max 


M [M ] 


1.400 


1.964 


1.400 


2.831 


2.822 


M B [M q ] 


1.778 


2.625 


1.767 


3.751 


3.756 


Q, [rads -1 ] 








4028. 


9547. 


9916. 


P [ms] 


_ 


_ 


1.560 


0.6581 


0.6337 


n c [fm -3 ] 


0.4842 


1.110 


0.4615 


0.7486 


0.8698 


p c [10 18 kgm" 3 ] 


0.7524 


2.059 


0.7125 


1.261 


1.517 


H c 


0.1747 


0.4514 


0.1588 


0.320 


0.370 


R cq [km] 


10.77 


10.71 


11.19 


16.54 


16.02 


r eq [km] 


8.573 


7.532 


8.946 


11.37 


10.86 


Tpolef ?~eq 


1. 


1. 


0.8935 


0.4618 


0.4762 


I [10 38 kgm 2 ] 


1.492 


2.155 


1.603 


6.534 


6.101 


J [GM|/c] 








0.7332 


7.084 


6.871 


cJ/(GM 2 ) 
Q/(MJ? 2 q ) 
c A Q/{G 2 M 3 ) 








0.3741 


0.8837 


0.8625 








0.0250 


0.0749 


0.0713 








0.7305 


1.172 


1.054 


T/W 








0.0348 


0.2100 


0.2007 


t/cq [C] 








0.1565 


0.5731 


0.5812 


z f 

-cq 


0.2742 


0.4768 


0.0463 


-0.3608 


-0.3634 


2 b 
'''eq 


0.2742 


0.4768 


0.5351 


2.584 


2.722 


■^polc 


0.2742 


0.4768 


0.2825 


0.8070 


0.847 


iVc 


0.6590 


0.4312 


0.6653 


0.4019 


0.374 










0.4203 


0.7320 


0.7540 


A c = B c 


1.439 


1.914 


1.432 


2.135 


2.234 


GRV2 


i. w-™ 


9. lO" 1 ^ 


1. 10" 4 


2. 10" 4 


9. 10" 4 


GRV3 


2. 10~ 12 


1. 10" 11 


2. 10 -4 


i. nr 4 


2. 10~ 3 



nally GRV2 and GRV 3 are the two virial error indicators 
discussed in Sect. 3.5. 



4-2. Evolutionary sequences 

A strange star which slowly loses energy and angular 
momentum via electromagnetic or gravitational radiation 
keeps its total baryon number constant. Therefore, we 
can compute evolutionary sequences of strange stars as 
sequences at fixed baryon mass Mb- Similarly to neu- 
tron stars, two categories of strange stars can be distin- 
guished: normal stars, which have a baryon mass lower 
than the maximum baryon mass of static configurations 



]\ ,j stat 
ivJ B,max 



= 2.625 B~ 1/2 M Q (cf. Table |l]), and 



supramas- 



sive stars, which have a baryon mass greater than M B ^* ax . 
Any normal star belongs to an evolutionary sequence 
which terminates at a static configuration. On the con- 
trary, supramassive stars exist only by virtue of rotation. 
The two families clearly appear in Fig. [|, which is a plot 
of the angular momentum as a function of the rotation 
frequency for evolutionary sequences. The supramassive 
sequences are not connected with the static limit 57 = 0. 
At both ends, they terminate by Keplerian configuration. 
The limiting case M B = M§ ^ ax is denoted by the thick 
line in Fig. ||. 

A characteristic feature of some supramassive se- 
quences also appears in Fig. |[ they can be spun-up by 



angular momentum loss: dil/dJ < 0. This effect is well 
known for Newtonian stars with soft EOS (adiabatic index 
close to 4/3) (Shapiro et al. 1990), as well as for (relativis- 
tic) neutron stars (Cook et al. 1994, Salgado et al 1994a). 
As can be seen from Fig. ||, this effect is very prononced 
for strange stars. 

Beside the fact that they can represent evolution path 
of rotating strange stars, another motivation for comput- 
ing constant baryon number configurations is that this 
permits a stability analysis, as discussed in Sect. 5.1 be- 
low. 

4-3. Comparison with previous works 

First of all, in the non-rotating case, we recover the re- 
sults presented in Haensel et al. (1986) for massless strange 
quarks [compare the M^f* model in Table [l] with Eq. (28) 
of Haensel et al. (1986)]. 

The first rapidly rotating model of a strange star 
has been computed by Friedman (unpublished, quoted in 
Glendenning (1989a, b)). The obtained maximum angular 
velocity is 6% lower than ours (last column in Table [l])^]. 
The corresponding gravitational mass is 7% lower than 



3 Friedman's result is presented for B 1/A := (hc) 3/4 B 1/4 = 
170 MeV, which corresponds to B — 108.7 MeVfm" 3 . It must 
be rescaled to B = 60 MeV fm~ 3 , according to the law (pp|), in 
order to be compared with results listed in Table [j]. 
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Fig. 5. Rotating strange star configurations in the mass - 
central density plane. The static limit and Keplerian limit 
are denoted by thick solid lines. Thin solid lines denotes 
sequences of constant baryon number; they are labelled 
by the baryon mass (in M Q units). The dashed line marks 
the limit of stability with respect to axisymmetric pertu- 
bations; configuration to the right of the dashed line are 
(secularly) unstable. Its intersection with the Keplerian 
limit curve is marked by an open square. It is the loca- 
tion of the fastest rotating stable star. Maximum mass 
configurations are denoted by circles 
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Fig. 6. Angular momentum as a function of the rotation 
frequency / — J7/(27r) along sequences of constant baryon 
number. Each sequence is labelled by its baryon mass in 
Mq units. The heavy black circles denotes Keplerian con- 
figurations. 



ours, and the central energy density is 7% greater than 
ours. 

The only other rapidly rotating model of strange star 
published in the literature is that of Lattimer et al. (1990). 
The maximum angular velocity obtained by these authors 
is 2% lower than ours (last column in Table 0) The 
corresponding mass is 12% lower than ours, the circum- 
ferential radius R cq 14% lower, the ratio T/W 10% lower, 
and the polar redshift (linked to the injection energy (3 
given by Lattimer et al. (1990) by z polc = - 1) 15% 

lower. The largest discrepancy arises for the baryon mass 
Mb- the value reported by Lattimer et al. (1990) is 21% 
lower than ours. However we suspect some error in the 
presentation of the results by Lattimer et al. (1990), be- 
cause (i) their ratio M-q/M is only 1.18, whereas it is as 
large as 1.34 for the non-rotating maximum mass model 
(for comparison, our value of Mb/M for the f2 max model 
is 1.35), (ii) in the neutron star case, large discrepancies 
with baryon masses presented in Lattimer et al. (1990) 
have been noticed previously by Salgado et al. (1994a) 
[see Sect. 5.1.2 of Salgado et al. (1994a)]. Note that the nu- 
merical code employed by Salgado et al. (1994a) has been 
successfully compared with codes from two other groups 
in Nozawa et al. (1998) and did not show any trouble with 
the baryon mass. Unfortunately, the baryon mass of Fried- 
man's model is not reported in Glendenning (1989a, b). 

Apart from the baryon mass problem discussed above, 
the discrepancy with Friedman and Lattimer et al. (1990) 
results remains rather large. It may be explained as fol- 
lows. Both Friedman's code [presented in Sect. II. a of 
Friedman et al. (1986)] and Lattimer et al. (1990) code em- 
ploy the Butterworth & Ipser (1976) finite difference tech- 
nique for numerically solving the partial differential equa- 
tions resulting from Einstein's equations. Butterworth & 
Ipser (1976) considered explicitly the case of a strong dis- 
continuity in the density field at the stellar surface [their 
Sect. III.c)ii)]. However they treat it only along the radial 
direction (by using a modified Lagrangian-polynomial fit 
in r) and keep a Legendre expansion for the angular vari- 
able /i := cos 9. This is correct as long as the stellar sur- 
face is exactly spherical, i.e. the star is static. But when 
the surface deviates from spherical symmetry (as is the 
case for rapidly rotating stars, cf. Figs. ^ and ||), there 
is also a discontinuity in the 9 direction. This disconti- 
nuity, when described by means of Legendre polynomials, 
which are smooth functions, inevitably generates spurious 
oscillations — the so-called Gibbs phenomenon. We be- 
lieve that it is this Gibbs phenomenon which explains the 
discrepancy between Friedman and Lattimer et al. (1990) 
results and ours. 



4 Lattimer et al. (1990) results (their Table 7) is presented 
for B = 1.634 x 56 = 91.5 MeVfrn" 3 . It must be rescaled to 
B = 60 MeVfm -3 , according to the law @, in order to be 
compared with results listed in Table 
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5. Maximally rotating configuration and scaling 
with B 

5.1. Maximally rotating stable configuration 

A natural bound on Q results from the mass shedding 
condition, corresponding to the Keplerian angular veloc- 



ity, f2x discussed in Sect. 3.4: stationary configurations of 
uniformly rotating strange stars with f2 > f2x do not exist. 
We will also require, that the stationary, uniformly rotat- 
ing configuration be stable with respect to axisymmetric 
perturbations (see e.g. Sect. 2.5.1 of Friedman 1998). We 
assume that the strange star is sufficiently hot, so that it 
is not subject to the r-mode instabilities (Madsen 1998). 
The maximum value of f2 for stationary, uniformly rotat- 
ing models of strange stars which satisfy these two con- 
ditions, will be denoted by f2 max , and the corresponding 
configuration will be called a maximally rotating one. 

The stability with respect to axisymmetric perturba- 
tions can be investigated by means of the turning point 
theorem established by Friedman et al. (1988). We use 
the following variant of it [also used by Baumgarte et al. 
(1998)]: along an evolutionary sequence (i.e. a sequence 
at constant baryon number, Sect. 4.2), the change in sta- 



bility is reached when the gravitational mass reaches a 
minimum. Using this criterion, we have determined the 
boundary separating stable and unstable configurations 
by linking the points of mimimum M along the constant 
baryon mass sequences in Fig. ||: the configurations lo- 
cated to the left (resp. to the right) of the dashed line in 
Fig. H are stable (resp. unstable) with respect to axisym- 
metric perturbations. Moreover, the intersection of this 
line with the fl = fix line in Fig. marks the maximum 
angular velocity allowable for stable rotating strange stars, 
i.e. the r2 max configuration. It is denoted by a square and 
is listed in the last column of Table [|. 

As noted by Cook et al. (1994) and Stergioulas & 
Friedman (1995), the f2 max configuration and the config- 
uration with maximum gravitational mass, M™* x , do not 
coincide, although they are close to each other for typi- 
cal neutron star models. We also found this for rotating 
strange stars, and the effect is more pronounced than in 
the case of neutron stars: it can be seen easily in Fig. 0. 
The parameters of the M™* x and r2 max configurations are 
listed in the last two columns in Table |l|. fi max is 4% higher 
than 0(M™* X ). Note that the M™* x configuration is on 
the stable side, so that the situation is similar to that 
presented in Fig. 9 of Stergioulas & Friedman (1995), or 
Fig. 8 of Koranda et al. (1997). 

5. 2. Formulas for f2 max 

The scaling properties of the equation of the stationary 
motion imply that fi max oc B 1 / 2 . The calculation of f2 max 
involves locating the threshold for the instability with 
respect to the axisymmetric perturbations, and requires 
therefore a higly precise numerical code for the 2-D cal- 



culations of stationary configurations. Therefore, the pre- 
cision with which the scaling formula for f2 max is fulfilled 
reflects the overall precision of our numerical calculations. 
Our numerical results for f2 max can summarized in exact 
formulae 

fl max = 9.92 10 3 rads" 1 , P min = ms . (24) 



We checked that the scaling il max oc v B holds with a very 
high precision. 

Within our simple model of EOS of strange matter, 
the constraints on the bag constant B, Sect. 2, imply 



0.513 ms < P min < 0.640 ms 



(25) 



In the case of dense baryon matter, exact results for the 
maximum rotation frequency of uniformly rotating neu- 
tron star models, for a very broad set of realistic causal 
EOS, show an interesting correlation with values of the 
mass and radius of the static configuration with maxi- 
mum allowable mass, M^*, Rff 1 (Haensel & Zdunik 
1989, Shapiro et al. 1989, Friedman et al. 1989, Friedman 
1989, Haensel et al. 1995). The most recent form of such 
an "empirical formula" for f2 max , derived in (Haensel et 
al. 1995), reads 



c- 



NS 



AP 



M 







10 km 



NS 



7730 rads" 



(26) 



where the prefactor CnSi which does not depend on the 
EOS of dense baryon matter, has been obtained via fitting 
"empirical formula" to exact numerical results for realis- 
tic, causal EOS of baryon matter. The "empirical formula" 
reproduces exact results with a surprisingly high precision, 
the relative deviations from Eq. (^6|) not exceeding 5%. In 
the case of strange stars, built of strange matter of mass- 
less quarks, the formula of the type ( |26| ) is exact, albeit 
with a different numerical prefactor. Using our numerical 
results, we get C ss = 7.84 10 3 so that C ss ^ C NS 
within better than 2%. Therefore, "empirical formula", 
derived originally for (ordinary) neutron stars, holds also, 
with very good precision, for strange stars. Our conclu- 
sion agrees with unpublished result of Friedman (quoted 
in Glendenning 1989a), and with Lattimer et al. (1990). 

5.3. Formula? for maximum mass configurations 

Rotation increases maximum allowable mass of strange 
stars, and the equatorial radius of the maximum mass con- 
figuration. Our results for the maximum mass of rotating 
configurations, M™^ x , and its equatorial radius, R eq! M^ x , 
can be summarized in two exact formulae, 



prot 



2.831 



Mr. 



km 



(27) 
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The ratio of the maximum mass of rotating configurations 
to that of static ones, and the ratio of the corresponding 
equatorial radii, are thus independent of B, 



AT 



max 



1.44 



prot 
ostat 



1.54 



(28) 



It has been pointed out by Lasota et al. (1996), that for 
realistic causal baryonic EOS the ratios M™* x /M^^ and 
i?™* Mmax /i?^ ax are very weakly dependent on the EOS, 
and can be very well approximated (within a few percent) 
by constants Cm = 1.18 and Cr = 1.34, respectively. 

In the case of strange stars the values of Cm and Cr 
are significantly higher; rotation increases the value of 
M max by 44% compared to 18% for neutron stars, while 
the corresponding increase in the equatorial radius is 54% 
compared to 34% for neutron stars. 

We conclude that "empirical formula" for f£ max has 
more universal character than the "empirical formulae" 
for M™^ x , -Req!M max ' proposed in Lasota et al. (1996); the 
former applies both for neutron and strange stars, while 
the latter describe neutron stars only. 



6. About triaxial instabilities of rapidly rotating 
strange stars 

6.1. Known results on the threshold in T/W for triaxial 
instabilities 

Rapidly rotating neutron stars can be subject to a sponta- 
neous breaking of their symmetry about the rotation axis, 
resulting from some triaxial instability, when their rota- 
tion rate exceeds a certain threshold (see e.g. Bonazzola 
& Gourgoulhon (1997), Friedman (1998) or Stergioulas 
(1998) for a review). In order for the symmetry breaking 
to take place, a (secular) dissipative mechanism has to be 
working. Basically two such mechanisms can be contem- 
plated: (i) viscosity and (ii) coupling with gravitational 
radiation (Chandrasekhar-Friedman-Schutz (CFS) insta- 
bility). For compact stars that are accelerated by accre- 
tion in X-ray binary systems, these triaxial instabilities 
may be an important source of gravitational waves, in 
the frequency band of the LIGO /VIRGO interferometric 
detectors currently under construction. Another situation 
where these triaxial instabilities might develop is in com- 
pact stars newly formed after a stellar core gravitational 
collapse. 

In the present paper, we have restricted ourselves to 
stationary and axisymmetric rotating strange star models. 
Therefore, we could not compute the triaxial instability 
threshold of rotating strange stars. However, an indicator 
of the stability of a rotating self-gravitating body is its ki- 
netic to gravitational potential energy ratio T/W, which 
can be computed for any equilibrium configuration. For in- 
stance, it is well known that a homogeneous Newtonian ro- 
tating body (MacLaurin spheroid) becomes secularly un- 
stable with respect to triaxial I — to = 2 perturbations 



(bar mode) if T/W > 0.1375 (Jacobi/Dedekind bifurca- 
tion point in the MacLaurin sequence). For compressible 
bodies (still in the Newtonian regime), this threshold can 
be lowered to 0.1275 (Bonazzola et al. 1996). The thresh- 
old of the CFS instability decreases with the mode number 
to for polar modes: for a Newtonian polytrope of abia- 
batic index 7 = 2, one has T/W > 0.079 for a to = 3 
polar mode, T/W > 0.058 for to = 4 and T/W > 0.046 
for to = 5 (Stergioulas & Friedman 1998). Higher order 
modes are likely to be damped out by viscosity. For axial 
r-modes, the CFS instability is generic, i.e. it occurs at 
any rotation rate (Andersson 1998). Madsen (1998) has 
however argued that the r-mode instability is suppressed 
by the quark matter bulk viscosity in new born strange 
stars. 

In the relativistic ratio T/W can be defined ac- 

cording to Friedman et al. (1986) prescription. This quan- 
tity reduces to the usual kinetic to gravitational energy 
ratio at the Newtonian limit, whereas its physical inter- 
pretation in the relativistic case is not so clear (in partic- 
ular the value of T/W is coordinate dependent). However 
it can be used to measure the importance of rotational 
effects. 

For the viscosity-driven instability, general relativis- 
tic terms increase the threshold on T/W: it becomes as 
high as T/W > 0.26 for a compactification parameter 
M/R — 0.2, according to the post-Newtonian analysis of 
Shapiro & Zane (1998). This stabilizing effect of general 
relativity onto the viscosity-driven bar instability has also 
been found by Bonazzola et al. (1998a). 

On the contrary, general relativity decreases the 
threshold of the CFS instability, down to T/W > 0.065 for 
to = 2 polar mode (this mode being always stable in the 
Newtonian regime for 7 = 2), T/W > 0.046 for to = 3, 
T/W > 0.035 for to = 4, and T/W > 0.029 for to = 5, for 
a 7 = 2 polytrope [Stergioulas & Friedman 1998, see also 
Yoshida & Eriguchi (1997)]. 

6.2. T/W ratio of rotating strange stars 

Having the above results in mind, we have computed the 
T/W ratio according to the prescription of Friedman et 
al. (1986), for each of our rotating stange star models. In 
particular, Fig. [7] shows the value of T/W as a function 
of the rotation frequency along evolutionary sequences. It 
appears clearly that the values of T/W for strange stars 
are much greater than for neutron stars (compare e.g. with 
Fig. 4 of Cook et al. (1994)). In particular, for the f2 ma x 
configuration, 



— = 0.201 , 



(29) 



(it does not depend on B) whereas it ranges from 0.105 to 
0.139 for the set of neutron star EOS examined in Nozawa 
et al. (1998). The fact that T/W for maximally rotating 
strange stars is significantly larger than that for neutron 
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Fig. 7. Kinetic to gravitational energy ratio T/W as a 
function of the rotation frequency / = f2/(27r) along se- 
quences of constant baryon number. Each sequence is la- 
belled by its baryon mass in Mq units. The heavy black 
circles denotes Keplerian configurations. 

stars, has been noted by Lattimer et al. (1990) and Colpi 
& Miller (1992). 

Moreover, contrary to neutron stars, the value of T/W 
remains very high, of the order of 0.2, for rapidly rotating 
strange stars of moderate mass: the value T/W — 0.2 
is reached at a rotation period P = 0.86 ms (resp. P = 
0.85 ms) for a baryon mass of Mb = 1.60 M Q (resp. Mb = 
1.80 M Q ), which corresponds to a gravitational mass M = 
1.31 Mq (resp. M = 1.47 M s ) (notice that in the static 
case Mb = 1.6, 1.8 M Q for normal sequences in Fig. [7] 
correspond to gravitational masses M = 1.27, 1.42 Mq, 
respectively) . 

All this seems to indicate that triaxial instabilities 
could develop more easily in rotating strange stars. How- 
ever, it has to be stressed that no actual stability analysis 
has been performed yet for strange stars which, in contrast 
to ordinary neutron stars, are bound not only by gravity, 
but also by confinement forces. These confinement forces 
are represented represented by the bag pressure, —B, act- 
ing on the strange star surface. The density profiles, Fig. 0, 
and adiabatic index 7, Fig. |l], are qualitatively different 
from those characteristic of rotating neutron star models, 
for which triaxial instabilities have been studied. 

7. Discussion and conclusions 

The numerical method used in the present paper is partic- 
ularly suitable for exact calculations of models of rapidly 
rotating strange stars. In contrast to numerical approaches 
used in previous calculations, the multi-domain method 
enabled us to treat exactly the huge density discontinu- 
ity at the surface of strange stars. We used three-domain 



grid in both r and 6 variables, adapted (adjusted) via a 
self-consistent procedure to the stationary shape of rotat- 
ing strange star surface. The precision of our calculations 
was checked using several tests. We found significant dif- 
ferences between our results, and those obtained by other 
authors in previous exact calculations of the models of 
maximally rotating strange star models. 

As noted in the early papers on the properties of 
strange stars, models of static massive (M <; M Q ) strange 
stars are characterized by global stellar parameters (mass, 
radius, moment of inertia, surface redshift) which are 
rather similar to those of ordinary neutron stars (Haensel 
et al. 1986, Alcock et al. 1986). In contrast to the static 
case, rapidly rotating strange stars exhibit qualitative dif- 
ferences with respect to rapidly rotating neutron stars. 

Strange stars rotating close to the break up (Keple- 
rian) frequency have kinetic to gravitational energy ratio 
T/W ~ 0.2, which is nearly twice larger than correspond- 
ing values for neutron stars. Moreover, this unusually high 
T/W ~ 0.2 is characteristic not exclusively for supermas- 
sive rotating models, but is typical also for normal rotating 
configurations with baryon mass lower than the maximum 
allowed mass for static models. In particular, T/W ~ 0.2 
is reached for the 1.4 M Q models rotating at P ~ 0.85 ms. 
Large value of T/W stems from a flat density profile com- 
bined with strong equatorial flattening of rapidly rotating 
strange stars. These particular features of rapidly rotat- 
ing strange stars are consequences of a special character of 
EOS of strange quark matter in the strange star interior, 
reflected by a specific density dependence of the adiabatic 
index. Strange matter inside massive (M > Mq) strange 
stars is very stiff in the outer (surface) layers and soft in 
the central core. Large value of T/W may signal a rel- 
ative softness (susceptibility) of rapidly rotating strange 
stars to triaxial instabilities. 

For simple strange matter EOS, used in the present 
paper, several exact scaling relations hold for extremal 
(maximum frequency, maximum mass) configurations of 
rotating strange stars. Maximum rotation frequency as 
well as maximum mass of rotating configurations exhibit 
an exact scaling with respect to the value of the bag con- 
stant. Empirical formula for f2 max for neutron stars, which 
relates this quantity to the mass and radius of static con- 
figuration with maximum allowable mass, holds with a 
high precision also for strange stars. 

Minimum period of rotating strange stars, stable with 
respect to axisymmetric perturbations, obtained for an 
acceptable range of the bag constant, is ~ 0.5 ms, and 
therefore similar to that characteristic of realistic mod- 
els of neutron stars. However, the effect of rapid rotation 
on the maximum allowable mass is significantly larger for 
strange stars than for neutron stars. Uniform rotation can 
increase the massimum allowable mass of strange stars by 
more than 40%, to be compared with only 20% increase 
for realistic models of neutron stars. 
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It should be stressed, that our numerical results for 
rapidly rotating strange stars, summarized in Eqs. (24,25), 
and Eqs. (27,28,29), have been obtained for a schematic 
EOS, Eq.(3). This particular EOS has been selected for 
the sake of simplicity and, last but not the least, for nu- 
merical convenience. The most general - generic - feature 
of an EOS of strange matter is the existence of a self- 
bound state at zero pressure, with Eq < 930.4 MeV. How- 
ever, the actual EOS of strange matter - should it exist - is 
expected to be different from that given by Eq. (3) . There- 
fore, results of the present paper neither represent the full 
range of possibilities for "theoretical" strange stars, nor 
they are expected to correspond to precise actual values 
of parameters of rapidly rotating strange stars - should 
they exist in Nature. 

The equation of state of strange matter, used in the 
present paper, has been derived neglecting strange quark 
mass and representing all effects of QCD interactions by 
one parameter - the bag constant B. However, this sim- 
plest (minimal) model exhibits generic properties of the 
EOS of a self-bound quark matter. Generally, inclusion of 
perturbative corrections, expressed in terms of the QCD 
coupling constant a c , and of the strange quark mass m s , 
imply shifting down and simultaneous narrowing of the 
window of the values of B compatible with strange mat- 
ter hypothesis (Farhi & Jaffe 1984). Typical values of m s 
are ~ 200 MeV/c 2 ; the effect of a c <> 0.3 on the EOS is 
then much smaller than that of m s . Scaling relations be- 
tween extremal strange star configurations, corresponding 
to different values of B (at fixed values of a c , m s ) are 
than no longer exact. In the case of static strange star 
models, they are still very precise (after the prefactors in 
these relations have been recalculated to include the effect 
of a c and m s ) and therefore useful (Haensel et al. 1986); 
similar situation is expected to take place in the case of 
rotating strange star models. In the static case, the effect 
of non-zero m s and a c consisted in making the maximum 
mass configurations less massive but more compact, and 
we expect similar effect for the extremal configurations of 
rotating strange stars. 

In this paper we restricted ourselves to the case of 
strange stars with superdense quark surface, built ex- 
clusively of strange quark matter (bare strange stars). 
In principle, one should consider also models of strange 
stars, covered by an envelope (crust) consisting of nu- 
clei immersed in electron gas (Alcock et al. 1986). The 
nuclei, forming a crystal lattice of the crust, are sepa- 
rated from strange matter by a repulsive coulomb barrier. 
Both coulomb barrier, and characteristic spatial gap be- 
tween nuclei and strange matter, result from an electric 
dipole layer on the quark surface, due to the nonuniform 
density distribution of electrons near the quark surface. 
The presence of electrons in strange matter, necessary for 
equilibrium of quark core core with crust, results from 
nonzero strange quark mass. The maximum mass of the 
crust on a M > M Q strange star was originally estimated 



as ~ 10~ 5 M (Alcock et al. 1986). Further studies of 
the crust-strange matter coexistence conditions led to an 
even lower value of the maximum mass of the crust on a 
M > Mq strange star, of the order of 10~ 6 M Q (Huang 
k Lu 1997). 

In the case of the static strange star models, the effect 
of the presence of the crust on the maximum allowable 
mass configuration is very small. Rapid rotation will in- 
crease the maximum mass of a strange star crust (Glen- 
denning & Weber 1992), but still its effect on the maxi- 
mally rotating and maximum mass configurations may be 
expected to be small: both f2 max and M™* x will be slightly 
reduced, as compared with the case of bare strange stars. 
It should be stressed, that the only existing calculation 
of rapid rotation of strange stars with crust (Glendcnning 
& Weber 1992) was performed within a version of slow 
rotation approximation. The multi-domain method with 
adaptable grid, used in the present paper, is also suitable 
for an exact treatment of rapid rotation of strange stars 
with crust, in which matter density exhibits a huge discon- 
tinuity at the quark matter-crust interface, with outward 
density drop by more than three orders of magnitude. This 
problem is now being investigated. 
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